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We study the Landau gauge gluon and ghost propagators of SU (3) gauge theory, employing the 
logarithmic definition for the lattice gluon fields and implementing the corresponding form of the 
Faddeev-Popov matrix. This is necessary in order to consistently compare lattice data for the bare 
propagators with that of higher-loop numerical stochastic perturbation theory (NSPT). In this paper 
we provide such a comparison, and introduce what is needed for an efficient lattice study. When 
comparing our data for the logarithmic definition to that of the standard lattice Landau gauge we 
clearly see the propagators to be multiplicatively related. The data of the associated ghost-gluon 
coupling matches up almost completely. For the explored lattice spacings and sizes discretization 
artifacts, finite-size and Gribov-copy effects are small. At weak coupling and large momentum, the 
bare propagators and the ghost-gluon coupling are seen to be approached by those of higher-order 
NSPT. 
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I. INTRODUCTION 



During the last years lattice studies on the Landau 
gauge gluon and ghost propagators have clearly not 
confirmed the asymptotic infrared behavior as postulated 
by the Gribov-Zwanziger scenario [a 6| or required by the 
Kugo-Ojima confinement criterion [7[. The inconsistency 
seems to be related to the treatment of the Gribov prob- 
lem, i.e., the non- uniqueness of the Landau gauge condi- 
tion in the non-pcrturbative regime (at least in finite vol- 
umes). In our opinion, it is still unclear, whether the lat- 
tice formalism — for a fixed gauge — requires a reformula- 
tion in order to become consistent with a BRST-invariant 
continuum theory or whether the continuum approach it- 
self needs to be reformulated (see also [1, Q ) . We would 
like to stress that until now no complete gauge-fixing pre- 
scription for Landau gauge has been agreed upon (see, 
e.g., [13, [HI and references therein). At the same time 
it is clear that the gauge-invariant lattice formulation of 
SU{N) Yang-Mills theories does provide confinement of 
quarks and gluons. Also in the context of gauge- variant 
approaches, there are clear signals that quarks and glu- 
ons (and ghosts) are not part of the physical spectrum, 
for example, by the violation of reflection positivity [T2I - 
[T3 |. The exact mechanism of confinement, however, is 
still unknown. It should be added though that a di- 
rect link between quark confinement, as measured by the 
Polyakov-loop order parameter, and the infrared behav- 
ior of ghost and gluon Green's functions has been estab- 
lished recently jl5j . 

However, the infrared asymptotics can only be one as- 
pect. Another, by no means less important, should be the 



computation of QCD's elementary two- and three point 
functions in the intermediate (around 1 GeV) and ultra- 
violet momentum region. This would allow for exam- 
ple to determine essential phenomenological parameters, 
like the QCD scale Ai^fg, or gluon and quark condensates 
(see, e.g., [l6l - l2]| ). Moreover, such calculations are im- 
portant to arrive at renormalized propagators, and even- 
tually also at vertex functions, which, calculated on the 
lattice and extrapolated to the continuum limit, could 
serve as input to a Bethe-Salpeter or Faddeev equations 
based hadron phenomenology (see, e.g., [22| for a status 
report). 

Therefore, the validity of multiplicative renormalizabil- 
ity in the nonperturbative regime as well as the speed of 
convergence to the continuum limit, and its uniqueness, 
are essential questions that need to be addressed on the 
lattice in this context. For the gluon field this has been 
done in the past [23l - l27| , in particular there with respect 
to the question of universality of its definition on the lat- 
tice. Also, more recently, first steps towards continuum- 
limit-extrapolated lattice data for the gluon and ghost 
propagators have been presented [1, [2^ . 

With our study we intend to provide further input to 
such projects, placing here particular emphasis on con- 
necting lattice Monte Carlo studies with those of lattice 
perturbation theory (LPT). Specifically we will confront 
its numerical variant, the numerical stochastic perturba- 
tion theory (NSPT) [H H^l, with lattice Monte Carlo 
(MC) data of the gluon and ghost propagators and the 
associated ghost-gluon coupling in Landau gauge [sij . 
This will help to quantify the region where predictions of 
NSPT are valid. 

So far, most MC studies of these gauge- variant objects 
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have used (what wc call below) the linear definition 



^(lin) 



traceless 



(1) 



with a, go denoting the lattice spacing and the bare cou- 
pling, respectively. On the lattice, however, the definition 
of the gluon field (and also of the Faddeev-Popov opera- 
tor) is in no way unique. One can of course equally well 
use another definitioUjfor example, that of the modified 
lattice Landau gauge |32|.^ or the quadratic definition 
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which has been studied and compared with the linear 
definition in Refs. [I^IIJ]. It was also shown there both 
converge towards the same continuum limit. 

Yet another definition is known as the logarithmic def- 
inition of the lattice gluon field. 
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This definition has been put forward for lattice MC stud- 
ies by Furui and Nakajima (see, e.g., [l^ HB-H^)- It is 
also the definition being used in lattice perturbation the- 
ory (see, e.g., the monograph [i^l) and NSPT. Since we 
are aiming at a quantitative comparison of NSPT with 
lattice MC data, below we will mainly concentrate on 
this definition. 

The respective control functional for the Landau 
gauge, the quadratic norm 

\\A\\'^J2 trA(x)At(x), (4) 
is given in terms of the lattice divergence 



(5) 

for each of these definitions. 

When expanded in terms of the lattice spacing all 
these definitions of the gluon field agree at leading order, 
but differ beyond that (see, e.g., [331). I* because of 
these differences why the corresponding Jacobian factors 
(Faddeev-Popov determinants) differ in the integration 
measure and why the respective lattice Landau gauge 



^ Note that the modified lattice Landau gauge (MLG) actually 
provides more than just another lattice discretization of the in- 
volved fields [33l . It is a potential candidate to overcome the 
0/0 Neuberger problem on the lattice |33|. |34| . Nonetheless, one 
could borrow the MLG lattice definitions for the gluon field and 
Faddeev-Popov operator if an alternative discretization of the 
standard lattice Landau gauge is desired. So far this would be 
possible for a U{1) or SU(2) gauge theory, and has been done 
for the latter for example in |35|| ). 



condition. (|| A| P) cannot be simultaneously satisfied (see, 
e.g., m, [231 or below). For a consistent setup of Landau 
gauge on the lattice one therefore has to employ the cor- 
responding lattice expressions for the gauge functional 
and the Faddeev-Popov operator. These are available 
in the literature for all above-mentioned approaches. For 
practical reasons, though, mostly the linear definition has 
been adopted. 

However, in order to assess the genuine non- 
perturbative effects in the measured two- and three-point 
functions, it is desirable to have as a reference point an 
understanding of the perturbative behavior of these func- 
tions in higher-order lattice perturbation theory. In re- 
cent years such higher-loop results for the lattice gluon 
and ghost propagators became available using iiumerical 
stochastic perturbation theory (NSPT) [2^|33|. These 
results are for individual momenta at a fixed lattice size 
and are usually obtained for the logarithmic definition. 
The advantage of NSPT is that the loop-order that can be 
achieved only depends on the available computational re- 
sources. There are no other restrictions. However, if one 
wants to confront the unrcnormalized NSPT results di- 
rectly with corresponding data from lattice Monte-Carlo 
(MC) simulations, one cannot expect the convergence of 
the former to the bare MC data, when different lattice 
gluon field definitions are employed. Of course, if mul- 
tiplicative renormalization is in place the bare propaga- 
tors for the different definitions are related through finite 
renormalizations, but for a comparison with NSPT this 
would bring additional uncertainties. For a quantitative 
comparison it is thus desirable to use the same lattice 
definition of the respective gluon field and the Faddeev- 
Popov operator. In addition, lattice NSPT results take 
the influence of the hypercubic group into account and 
thus allow for a comparison of the respective propagators 
at more off-diagonal momenta which are usually excluded 
by cuts. 

Wc therefore find it worth to complement the existing 
data for the Landau gauge gluon and ghost propagators 
by new sets for the logarithmic definition of the lattice 
gluon fields and the corresponding Faddeev-Popov op- 
erator. The bare propagators can then be confronted 
directly (i.e., without any additional renormalization) to 
the results from NSPT. We think that this way much 
more precise information can be obtained on the regime 
where NSPT holds. 

Note that for a renormalization-group invariant, like 
for example for the ghost-gluon c oupling , denoted 
a^"{q'^) in the minimal MOM scheme [il|4l|], no rescal- 
ing of the data would be necessary, when comparing data 
obtained for different lattice discretization.^ Data for this 
coupling, extracted either for the linear or the logarith- 



It is important though, that the bare lattice gluon and ghost 
dressing functions is for the same type of lattice discretization, 
and the correct tree-level value is obtained when the links are set 
to one (see also ISSlI ). 
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mic definition, should match up completely, apart from 
artifacts due to the lattice discretization. This has been 
argued and, for SU{2), also exphcitly shown already in 
[35|- This coupling thus serves a reasonable object to 
systematically investigate discretization effects. Below 
we will further corroborate these findings for the gauge 
group SU{3) and compare the available data sets to cor- 
responding ones of NSPT. 

Altogether, our aim is threefold: 

1 . To check universality of the employed lattice defi- 
nitions by comparing MC results obtained for the 
logarithmic and linear definition. Such a study is 
not really new, but we do it simultaneously for 
the gluon and ghost propagator and the associated 
ghost-gluon coupling. We confirm that, at least in 
the momentum ranges considered throughout this 
paper, the propagators are indeed related to each 
other by multiplicative renormalization, and that 
these renormalization constants are such that they 
exactly cancel each other in the coupling as ex- 
pected. 

2. To compare nonpcrturbative and perturbative lat- 
tice results obtained for exactly the same set of pa- 
rameters. We emphasize that (in this study) we 
are not intending to fit lattice results in the pertur- 
bative range with higher-order perturbation the- 
ory as obtained in the continuum by [4]1 - |43j . This 
has been the intention in the lattice references 

[21I . . Here we are rather interested to locate the 
momentum region where predictions of NSPT are 
still a good approximation to Monte Carlo results. 

3. To report in some detail on the realization and per- 
formance of the gauge-fixing algorithm required for 
the logarithmic definition, as well as to assess the 
importance of the Gribov ambiguity in this setting. 

The paper is organized as follows: In Section |ll] we 
briefly review lattice Landau gauge for the linear and the 
logarithmic definition of gluon fields. Section Hill is then 
devoted to gauge-fixing algorithms, first for ^i'+^/2p' 
and then for A^^°^} ,„ . For the latter we test three 

x+fj,/2,fj. 

different implementations: an unaccelerated, a Fourier- 
accelerated and a multigrid-accelerated gauge-fixing al- 
gorithm. We will argue that one should preconditioning 
these by first fixing the gauge field configuration such 
that ^^'"^/2 fj. transversal, before the actual gauge- 
fixing of ^^'^^^/2 starts. A comparison of the three 
different implementations including their parameters is 
given at the end of Section Hill In Section HVl we present 
expressions for the gluon and ghost propagators for the 
logarithmic definition. A brief study on the impact of 
the Gribov ambiguity is also discussed there. Further 
results are then presented in Section |Vl where we first 
compare the gluon and ghost propagators for the log- 
arithmic definition with those for the linear definition, 



and then discuss lattice discretization and finite-volume 
effects. We will demonstrate that data for the coupling 
a^'^(g^) matches up for both definitions without any 
rescaling. In Section IVll we finally confront our MC re- 
sults for the propagators and the coupling to the recent 
resuhs from NSPT [HllOl. We will argue that for large 
/3 = 6/(7o and large momenta the MC data must be re- 
stricted to the trivial (i.e. real-valued) Polyakov loop sec- 
tor in order to reach good agreement.'^ In Section IVIII we 
will draw our conclusions. A detailed discussion of the 
multigrid-accelerated gradient algorithm follows in the 
Appendix. 

II. LATTICE IMPLEMENTATIONS OF 
LANDAU GAUGE 

On the lattice the fundamental degrees of freedom are 
the link variables Ux,fi; which are elements of the SU{3) 
gauge group. If one is interested in gauge- variant quan- 
tities, like for example the gluon propagator, one has to 
adopt a definition for the gluon field ^2;+/i/2,/i in terms of 
these links. Above, we mentioned three ways of defining 
a gluon field on the lattice, the linear and logarithmic 
definition will be considered below. 



A. Linear definition of the gluon field 

We first recall the linear definition [Eq. ([ij]. For this 
the Landau gauge condition on the lattice (setting the 
lattice spacing a = f ), 



A("")(x) = 




is realized, if the gauge functional 

^f'[5]-i^E (1-5 ^etr^O (7) 

is in a (local) minimum for a given gauge-fixed configu- 
ration 

'u,,, = g.u^^Mi^ ■ (8) 

Here, Qx G SU{3) is the gauge transformation field that - 
finally - should put the unfixed gauge field Ux.^ to Lan- 
dau gauge. As there are many local minima ("Gribov 
copies" ) we assume that unique gauge fixing is achieved 
by searching for the global minimum for each configu- 
ration U. In practice, one can only try to get as close 



3 Note this was argued for SU(2) already in liSl . 
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as possible to the global extremum. We call this general 
prescription "minimal Landau gauge" . The minimization 
of F^''"'' [g] can be accomplished by an overrelaxation al- 
gorithm or a combination of a simulated annealing and 
overrelaxation (see below). This allows to fulfil the Lan- 
dau gauge condition ([6]) with the required local numerical 
precision. 



B. Logarithmic definition of the gluon field 

In the continuum, Landau gauge can also be formu- 
lated as a minimization problem of the functional 

J^i™"*^[<?] = ]^ E / I'M^) '^m(^)] (9) 



with respect to g{x) £ SU{3) acting on Afj_{x) according 
to 

<^A^{x) = g{x)A^{x)g\x) + - g{x) d^g\x) . (10) 

90 

The direct transcription of the continuum extremization 
problem to the lattice leads to the minimization of a lat- 
tice gauge functional 
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for the logarithmic definition of the lattice gluon field 
[Eq. ([3])] . The violation of transversality then can always 
be checked by computing the divergence 



(log) 



-A 



(log) 



^0. (12) 



Note that the logarithmic definition requires a diago- 
nalization of the neighboring SU{i) link matrices each 
time the divergence is evaluated. 

We will see that the two lattice Landau gauge con- 
ditions [Eqs.® and (jl2p ] - if imposed on an arbitrary 
gauge field configuration - will result in rather different 
gauge-fixed fields. However, the gluon propagators cal- 
culated from the respective gluon fields, i.e, A^^ 



A'"^"^^. ,„ , will be seen to be related to each other ap- 
proximately by a finite multiplicative renormalization. 
Note also that the gluon fields transformed into momen- 
tum space will be transversal, q^Afj^{q) = 0, only if the 
gauge-fixing procedure fits to the respective definition of 



(log) 



A. 



III. LANDAU GAUGE FIXING: DIFFERENT 
ALGORITHMS ARE REQUIRED 

A. Linear definition 

The linear form of the gauge functional in Eq. ([7]) sug- 
gests to use a relaxation method for its minimization (Los 
Alamos type gauge-fixing): at given lattice site x the 
gauge transformation field gx is replaced by g'^ such that 
the expression 



$netr{5i Wx} 
is maximized for a given Wx, where 



(13) 



This is also known as a "projection onto 5L/(3)" 



g'x = Proj 



sum 



(14) 



For SU{3) g'x is typically found using the Cabibbo- 
Marinari decomposition [45[, such that the local SU{3) 
update of gx proceeds via three successive SU{2) up- 
dates. Formally, the update can be viewed as the re- 
placement 



9x-^ g'x = rx Qx with = g'^ gl 



(15) 



The speed of convergence is usually improved by replac- 
ing the relaxation steps by overrelaxation (OR) steps. 



9x -^r^ gx. 



(16) 



The overrelaxation parameter w has to be optimized in 
the interval 1 < w < 2. The required power of the update 
matrix Tx is approximated via the truncated series 
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E 



where = 3 or 4, and 



[rx ~ If , 



7n(w) = 



r(^ + i) 
r(w + 1 - n) 



(17) 



(18) 



In order to check the Gribov ambiguity, the OR pro- 
cedure can be repeated for a number of initial random 



gauges g^^lti^i 



, -^copy typically resulting in dif- 

,(') 



x+fi/2,fj.- 



ferent final gauge transformations gf^^^i for a given lattice 
gauge field U. The different gauge-transformed fields are 
known as "gauge copies" . They lead to a corresponding 
distribution of gauge-functional values F^*^ = -F^''"'' [g^*^] 
(the set of local minima) instead of a unique, but usually 
unknown absolute minimum. 

To shift the distribution of F^'' to smaller values, the 
initial random gauge transformation can be replaced by 
one obtained using a simulated annealing algorithm (SA) 
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for gauge-fixing. SA is a Markov chain MC process sim- 
ulating a Gibbs measure of the form 



Pu.tM oc cxp 



To 



(19) 



for the field of gauge transformations g, where the gauge 
temperature Tg is lowered step by step after one or a few 
Markov steps according to some protocol ( "SA schedule" ) 
such that the process never passes through real equilib- 
rium states. To our knowledge, the method has been very 
successfully applied for the first time fixing to maximally 
Abehan gauge in (46j . 

It turned out that with respect to computer time as 
well as with respect to finding smaller F^'^\ repeating 
a combined SA+OR algorithm is more efficient, than 
repeating the OR algorithm with initial random gauge 
transformation (see [43, 13 ) ■ For the results presented 
below we have employed the SA-I-OR algorithm to gauge- 
fixing the fields A^'™) [Eq. ([T])] but also to precondition 
the gauge-fixing of [Eq. dSj] (for details see below). 

For the SA algorithm we have changed the gauge tem- 
perature Tg{is) for every new MC sweep (is) according 
to the SA schedule 



(T, 
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iter 
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(20) 



proposed in 147[ with a restriction to Niter = 3500 itera- 
tions. We have used the initial maximal gauge temper- 



ature value To 



0.45 and the final minimal value 



0.01. 



B. Logarithmic definition 

For the logarithmic definition of the gluon field (for 
brevity we call the corresponding algorithm logarithmic 
gauge fixing), the starting point is the local divergence 
Eq. (fT2)) evaluated at all lattice sites x. The gauge trans- 
formations (at a given lattice site x) is updated locally 
by gx ~^ fxQx where is the exponentiated local diver- 
gence of ^A^ evaluated at x according to [Eq. P^ ]. i.e., 



exp -za (^a/4'°s))(a;) 



(21) 



instead. Here F denotes the Fourier transformation from 
the space-time lattice to the 4-momentum lattice and 
the reverse transformation, respectively (for see 
and ([5n)) V One recognizes a sort of smearing of the 
divergence with the inverse Laplacian. We notice that 
Furui and Nakajima [H^l have also used a non-local algo- 
rithm in the form of a Newton-Raphson-typc construc- 
tion of Tx using the non- vanishing divergence as a source 
and the inverse Hessian (Faddeev-Popov operator, ex- 
panded in to some finite order) in place of the inverse 
Laplacian. 

Fast Fourier transformations are less efficient if the lat- 
tice is fully parallelized. This could become a problem 
for large lattices. We circumvented this problem by us- 
ing the multigrid- accelerated steepest gradient algorithm 
5lj . In fact, the last expression for Tx can be written as 



cxp 



'Zmax 



A-i(^a^M('°s))(x) 



(23) 



The inversion of the Laplacian on an arbitrary source is 
a standard problem for a multigrid algorithm. Details on 
our implementation of the multigrid-accelerated gradient 
algorithm arc given in the Appendix. 

We faced the problem that the logarithmic gauge fixing 
fails to finish successfully for 10% to 50% of all attempts, 
when starting from a random gauge transformation. The 
exact failure rate depends on the gauge-fixing method, 
and grows with the lattice size. A similar problem only 
rarely occurs in case of the linear gauge fixing using the 
standard OR algorithm. We found, however, that the 
logarithmic algorithm works successfully in 99% of the 
cases if the divergence [Eq. (IT^ j is already sufficiently 
small at the start. This can be achieved by precondi- 
tioning the logarithmic by a linear gauge fixing. The 
remaining 1% of cases is dealt with by just repeating the 
gauge fixing. 



C. Comparing the two gauge fixing prescriptions 

For both types of gauge fixing a stopping criterion is 
needed. This criterion is fulfilled as soon as the respec- 
tive divergence, {d^^A^^^^''^ {x) or {x) is suf- 
ficiently small. We have applied the criterion 



The step size a has to be tuned. We call this method the 
local unaccelerated steepest gradient algorithm (the orig- 
inal "Cornell type" gauge fixing discussed in Ref. (49|). 

It is well-known that this update suffers from critical 
slowing down which can be ameliorated using a Fourier- 
accelerated version [49| '. at each lattice site gx — ^ i^xQx 
using 



cxp 
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(^a/A(;°s))(^) 



(22) 



max tr 



(9/A,)t [x) {d,<^A,) [x) 



< 10 



-14 



(24) 



To demonstrate the relation between the two different 
gauge-fixing algorithms we compare now the history of 
the overrelaxation and of the Fourier-accelerated gradient 
algorithm. Both are applied to the same configuration on 
a 16* lattice (generated using the Wilson plaquettc action 
at /3 = 6.0). 

We used w = 1.68 and = 3 [sec Eq. ^] for the 
linear gauge-fixing method (in this case the OR algo- 
rithm without the SA-prcconditioning step) and a = 0.07 
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FIG. 1: History of maximal values of the squared divergence of the two gluon field definitions during gauge fixing by overre- 
laxation (left) and the Fourier-accelerated gradient method (right). Example configuration on a 16"* lattice at /3 = 6.0. 
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FIG. 2: Same as in Fig. [T] but for the space-time average of the squared divergence of the respective Afj,. 



[Eq. (|22p ] for the logarithmic gauge fixing (i.e., the 
Fourier-accelerated gradient algorithm). 

In Fig. [1] we compare how the local maxima of the 
squared divergence of A^''") and ^('°s) behave while ei- 
ther the overrelaxation (left panel) or Fourier-accelerated 
gradient algorithm (right panel) proceeds. When the OR 
algorithm fixes the Landau-gauge the divergence of A*^''") 
is successively reduced. This cannot be simultaneously 
achieved by A'-^°^K On the other hand, when the gradi- 
ent algorithm fixes the Landau-gauge the divergence of 
^(log) becomes smaller with every iteration, while the di- 
vergence of A*^''"^ stays almost unchanged. We notice 
that the local maximum of the squared divergence of 
the "wrong gauge field" never becomes less than 0(1). 
This is the result of only few isolated defects, as it can 
be seen by a comparison with the corresponding space- 
time averaged quantities (Fig. [5]) . The space-time aver- 
age defined through ||A|p [Eq. (|1])] reaches a precision 
of 10"'^, typically after O(IOO) or 0(200) iterations. To 
reach the same precision for all lattice sites more itera- 
tions are needed. We notice that the space-time aver- 
age of the squared divergence of the "wrong gauge field" 
never becomes less than 0(10^^). Comparing the linear 
and the quadratic definition of the gluon field the au- 
thors of Refs. [1^ have observed a similar difference 
of ||A(<i^'^<i)||2 and ||A("")||2. 

Local operators constructed in terms of the gauge- 



fixed gluon fields can only be sufficiently precise if the 
respective gauge-fixing method has been used for each 
Afj_. However, as a preconditioner we can (and wc do) 
use the linear gauge-fixing (via the SA-I-OR algorithm) 
before applying logarithmic gauge fixing. As mentioned 
above, this helps much to reduce the number of unsuc- 
cessful gauge-fixing attempt for the latter. For the linear 
definition the SA+OR algorithm always converges. 

Preconditioning is not without effect on the final gauge 
fixing, however. The quality of gauge fixing achieved 
by the (linear) preconditioner infiuences the quality of 
the final (logarithmic) gauge fixing. This is illustrated 
in Figs. 13 and H at /? = 6.0 for a 16'' and 32^ lattice, 
respectively. These figures show the effect of replacing 
pure OR by SA-I-OR (with respect to A*^''"^) as precon- 
ditioner for the logarithmic gauge fixing (performed with 
the multigrid-accelerated gradient method) on the dis- 
tribution of the final gauge functional values F^'"^"* [5 '■'■']. 
This comparison is again for a. 16"^ lattice at /? — 6.0. 

As well-known for the linear gauge fixing, we typically 
find smaller gauge-functional values also for the logarith- 
mic case, when the SA-I-OR algorithm is used instead of 
the OR algorithm for this preconditioning step. This be- 
comes even more pronounced increasing the lattice size 
where the number of local minima naturally increases. 
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FIG. 3: The effect of replacing OR (left) by SA+OR (right) as preconditioner. The figure shows a histogram of final i^('°s) 
values for 100 gauge-fixed copies (16* lattice, /3 = 6.0). The multigrid-accelerated gradient method was used for the gauge-fixing. 
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FIG. 4: Same as Fig. El but for a 32* lattice. 



D. Performance of different realizations of the 
logarithmic gauge fixing 



We compare now the logarithmic gauge fixing (without 
SA+OR preconditioning) in three versions: the unaccel- 
cratcd gradient algorithm, the Fourier-accelerated gra- 
dient algorithm and the multigrid-accelerated gradient 
algorithm. The advantage of the latter is that it is easy 
parallelizable. Our code employes a further developed 
version of the algorithm used by Cucchieri and Mendes 
for Landau gauge in SU (2) gauge theory [52j. In the Ap- 
pendix our implementation of the multigrid-accelerated 
gradient algorithm is described in more detail. 

In Fig. [5] (left) we show the average number of itera- 
tions of the multigrid-accelerated algorithm as a function 
of the step size parameter a (again for a 16* lattices with 
/3 = 6.0). With increasing a the mean number of itera- 
tions decreases monotonously until it reaches an optimal 
value a = aopt- Further increasing a beyond that value 
leads to instabilities and is therefore avoided. In Table U 
we summarize our values on ofopt , for the three algorithms 



and different lattice sizes and /?. Due to limited comput- 
ing resources, for the bigger lattices we could afford to 
fix the gauge only by means of the multigrid-accelerated 



TABLE I; Optimal values aopt for the three different logarith- 
mic gauge-fixing algorithms for /3 = 6.0 and 9.0, and different 
lattice sizes. Due to the high demand of computation time, 
Qfopt has not been determined for the unaccelerated and the 
Fourier-accelerated algorithm on lattices larger than 16*. 



algorithm 


/? 


8* 


12* 


16* 


24* 32* 


unaccelerated 


6.0 


0.130 


0.130 


0.110 




Fourier-accelerated 


6.0 


0.070 


0.070 


0.070 




multigrid-accelerated 


6.0 


0.070 


0.070 


0.075 


0.70 0.70 


unaccelerated 


9.0 


0.125 


0.125 


0.125 




Fourier-accelerated 


9.0 


0.065 


0.065 


0.065 




multigrid-accelerated 


9.0 


0.065 


0.065 


0.075 
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FIG. 5; The number of iterations needed to reach the stopping criterion. Left: as a function of the step size parameter a of 
the multigrid-accelerated algorithm (lattice size le"*, /3 = 6.0). Right: as a function of the lattice size for the three investigated 
logarithmic gauge fixing algorithms (/3 = 6.0). 



TABLE II: The exponent and coefflcient, z and C, of Eq. ([25l 
for the three different logarithmic gauge-fixing algorithms. 



algorithm 


C 


z 


unaccelerated 


110 


0.33(4) 


Fourier-accelerated 


13 


0.3(4) 


multigrid-accelerated 


11 


0.3(4) 



gradient algorithm (in its parallelized version).^ Fig. [5] 
(right) presents the scaling of the average iteration num- 
ber with the lattice size for the three logarithmic gauge- 
fixing algorithms with their respective ctopt (again all for 
/3 = 6.0). We find that these numbers to scale like 



Mtcr ^CV' 



(25) 



Values for C and z are summarized in Table ITll 

Apparently, the two accelerated algorithms perform 
much better compared to the unaccelerated one, which is 
mainly due to a much smaller C [compare 0(10) versus 
0(100)]. The values for z arc almost the same for all the 
three algorithms. We favor the multigrid gauge-fixing 
algorithm because it is easy to parallelize. 



IV. GLUON AND GHOST PROPAGATORS FOR 

ylL""' AND a!,'"*^' 



A. Gluon propagator 

We are interested in the gluon and ghost propagators 
for the linear and logarithmic definition. With the lattice 
gluon field A'''") and ^ respectively, the bare gluon 



The lattice ensembles at /3 = 9.0 have only been used for com- 
parison with NSPT. Larger lattice sizes could be simulated by 
NSPT only quite recently. They were not available at the time 
when the MC studies described here were finished. 



propagator is defined as 



D 



ah 



(26) 



where (...)[; denotes the ensemble average over gauge- 
fixed configurations. As in most of the applications we 
evaluate it within the lattice Fourier representation 



with 



E 



e 



(27) 



(28) 



where the abbreviation A: • x — V], 



has been used. 



^ denotes the lattice size in /i-dircction. The inte- 
gers k^ G (— L^/2, +L,,/2] count the momentum modes 
within the Brillouin zone. The lattice momenta can be 
written in two forms. 



and 



(It^ik) 



2 Trfc 
— sm — ' 
a L, 



2 . a9p 
— sm — - 
a 2 



(29) 



(30) 



The latter is the one that appears in the lattice tree-level 
expression for the gluon propagator, and therefore taken 
as the corresponding physical momentum. 

For the lattice spacing dependence a(/3) we adopt [ssj 
and use ro = 0.5 fm to assign physical units to a. Ta- 
bic mil lists the lattice spacing values used in this study. 

Supposed all non-diagonal components in the color 
(a, b) and the Euclidean indices (/i, v) vanish, one can 



average over the diagonal elements 



(31) 
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TABLE III: Values for the lattice spacing a(/?) and its inverse 
as used in this study. We used the formula given in [s^] with 
ro = 0.5 fm. 





a(/3) in fm 


a-^{l3) in GeV 


5.8 


0.1364 


1.4464 


6.0 


0.0932 


2.1184 


6.2 


0.0677 


2.9137 


6.4 


0.0513 


3.8445 



where (written up to fourth order in the gluon field) 



1 - u, 



cab 



t 



2 



(37) 



\ / 2 \a'b 1/^ \ab 

12 I^^^+I^mJ " 720 [^^+%^^) 



The factor 1/3 is due to the transversahty of the gluon 
field with respect to the vector g^(fc) that leaves only 
three independent modes. The gluon dressing function is 



(32) 



B. Ghost propagator 

For the ghost propagator the situation is somewhat 
different. It is the two-point function of the ghost fields 
c° and c^, and these fields are only implicitly defined 
through the Faddeev-Popov operator M°^. Therefore, 
the bare ghost propagator is defined as the ensemble av- 
erage of the inverse Faddeev-Popov operator, i.e.. 



v/u 



(33) 



As mentioned the form of the Faddeev-Popov operator 
depends on the definition adopted for the gluon field. 
Generally, on the lattice one can decompose the Faddeev- 
Popov operator as follows: 



j^jab 



Af 5. 



xy 



Sx 



Cab 



Sx^ 



(34) 



For the linearly defined gluon fields one determines the 
Faddeev-Popov operator as the Hessian of the gauge 
functional with respect to infinitesimal gauge trans- 
formations. This leads to the following form of the ad- 
joint representation matrices entering p4p 



Af = 5Hetr 



{T-,T''}J2iUx„ 



+ Ux-f.^t.) 



2 9^etr[r" C/,_^,^] 



(35) 



a 



ah 



where (a — 1, ... ,S) denote the generators of SU{3). 

For the logarithmic definition the following form of the 
adjoint representation matrices is used, 



ab 



A 



r>ab 
Bx,,j. 



c: 



ab 

X,fJ. 



''''x-iJ.,11 



ab 



abc 



(36) 



A'= 



jabc 



with the gluon field A"'', its square (A^)"'' etc. taken in 
the adjoint representation. 

The Faddeev-Popov operator is inverted with a 
Laplacian-preconditioned conjugate gradient algorithm 
and a color-diagonal plane wave source as explained 
in m. 

For both the linear [Eq. ([55)1 ] and the logarithmic defi- 
nition [Eq. p7p] the ghost propagator in the momentum 
space is given by [q = q{k)\ 



(38) 



The corresponding ghost dressing function is 

ZGh{q^)^q^G{q^). (39) 

C. Momentum cuts 

On the lattice, the 0(4)-symmetry of the continuum 
Euclidean space-time is broken to the discrete H{A)- 
symmetry. In order to minimize the resulting artifacts, 
due to the finite lattice volume and spacing we app ly 
two momentum cuts: a cone and a cylinder cut |55| . ^ 
The cone cut removes all (low) lattice momenta k — 
(fci, ^2, fca, ^4) with at least one vanishing k^. It is applied 
to minimize finite-volume effects associated with these 
lattice momenta. The cylinder cut removes all momenta 
which are not close to a multiple of one of the space-time 
diagonal unit vectors n = (1/2)(±1, ±1, ±1, ±1). For a 
symmetric lattice, this criterion can be formulated as 



< 1 



(40) 



These two cuts removed most of the lattice artifacts from 
the data. The remaining artifacts can be assessed using 
different lattice spacings a and lattice volumes V . 



For an alternative approach see, e.g., Ref. [56 
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FIG. 6: Left: the bare gluon propagator D{q^) from the logarithmic gluon field definition and the relative difference SD{(f) 
in percent between fc- and bc-configurations generated at /3 = 6.0 for lattice size 32*. Right: the corresponding bare ghost 
propagator G{q^) with its relative difference 5G{q^). 



D. Gribov ambiguity for -propagators 



In order to carry out a first check of the Gribov am- 
biguity of the gluon and ghost propagators for the loga- 
rithmic definition aII"®-* we used 100 (30) configurations 
for a 16" (32^) lattice (again all for /3 = 6.0). For the 
larger lattice we expect a stronger influence of the actual 
selection among the gauge-fixed copies. 

Employing the SA-f-OR preconditioned, multigrid- 
accelcratcd gradient algorithm we fix the gauge for each 
configuration 10 times starting from different random ini- 
tializations of the gauge transformation field g^. As in 
[13 for each configuration we have calculated the gluon- 
and ghost propagator for the first (random) gauge-fixing 
attempt (first copy, fc) and for the gauge field with the 
lowest gauge functional value achieved (best copy, be). 
Fig. [5] shows the gluon propagator D{q^) and the ghost 
propagator G{q'^), respectively, averaged over the first 
and best copies, respectively. In the lower panels of the 
figures we present also the corresponding relative differ- 
ences SD and SG = "((fc— bc)/bc)" of fc- and bc-results 
given in percent. 

The gluon propagator docs not show any effect of 
the Gribov ambiguity beyond statistical noise ("Gribov 
noise") over the whole momentum range, while the ghost 
propagator seems to exhibit a slight systematic shift 
within the low momentum region. This shift is approxi- 
mately 1% for the 32^ lattice while for the 16^ lattice it 
turns out to be negligible. The small effect of the Gribov 
ambiguity is certainly a consequence of the precondition- 
ing step for which the SA-I-OR algorithm has been em- 
ployed. Whether it becomes more enhanced, when taking 
global Z(3)-flip transformations into account P. [ssl. [soj . 
remains to be seen. 

In the following we will always rely on bc-results. 



TABLE IV: Statistics of Monte Carlo ensembles. A'^conf gives 
the number of analyzed configurations and A'^copy of inspected 
gauge-copies for each. Values are the same for all /3. 



lattice 


Aconf 


copy 


12* 


200 


10 


16* 


100 


10 


24* 


50 


10 


32* 


30 


10 



MONTE CARLO RESULTS FOR THE 
PROPAGATORS 

A. Comparing propagator results for 



A 



(lin) 



and A 



(log) 



To get better insight into discretization effects we com- 
pare now the gluon and ghost dressing functions for the 
logarithmic and linear definition. Unfortunately though, 
we have to restrict the discussion to relatively small lat- 
tice sizes. 

The data presented arc based on ensembles of gauge 
field configurations with statistics as given in Table IIVI 

We consider first the gluon dressing function calculated 
for the linear and logarithmic definition on a and 16^ 
lattice at (3 — 6.0 and /3 = 9.0. Note that the latter was 
chosen only in order to compare with available NSPT re- 
sults (see Section |VI|. Fig.[7]shows the data for the bare 
dressing function versus the lattice momentum squared. 
We clearly see the expected momentum- independent off- 
set between the results for the logarithmic and the linear 
definition. 

This is even better seen in Fig. [5] where the ratio 



C'gluon('7^) 



D(lin)( 



q 



D(i°g)(g2) 



(41) 
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FIG. 7: The bare gluon dressing function for /3 = 6.0 (left) and /? — 9.0 (right) versus the lattice momentum squared. The 
data is for two lattice sizes. Filled symbols are for the logarithmic definition and open symbols for the linear. 



is shown versus momentum. For both /3-values, we ob- 
serve the ratio to be constant within statistical errors over 
the whole momentum region, and to depend on /3. Thus 
the bare gluon propagators differ for the two definitions 
but are related to each other by a finite /3-dependent mul- 
tiplicative renormalization constant. As a consequence 
both definitions will lead to the same propagator when 
renormalized in a MOM scheme, the latter being defined 
by the condition that the propagator equals its tree-level 
expression at some subtraction momentum ^ fi^. Of 
course, this multiplicative renormalizability has been nu- 
merically demonstrated here only for finite volume and 
corresponding restricted momentum range under consid- 
eration. 

Also for ghost propagator we clearly see the offset be- 
tween the dressing functions for the logarithmic and lin- 
ear definitions (see Fig. |9] for the data at /3 = 6.0 and 
9.0). 

Similar to Fig. |51 in Fig. [TO] we show the ratio 



Cghost('? ) 



2, G("")(g2) 



G(log)(g2 



(42) 



as a function of the momentum squared . For both val- 
ues of /3, we see an approximately constant ratio over a 
wide momentum range. The deviation seen at the small- 
est momenta for /3 = 9.0 remains within statistical errors. 
In Table |V] we list the values for Cgiuon and Cghost • As 
expected, their ratios happen to be related as 



^2 

Cghost 



1. 



(43) 



This implies that the ghost- gluon coupling a^^(g^) (see 
Section IV C|) determined directly from the gluon and 
ghost dressing functions is the same for the logarithmic 
or linear definition. 



TABLE V: The ratios Cgiuon and Cghost [Eqs. (gl]) and ^ 
for the linear and logarithmic definition. 



lattice 


P 


Cgiuon 


C'ghost 


12'* 


6.0 


0.82 ± 0.02 


1.1013 ±0.0007 


12* 


9.0 


0.91 ±0.03 


1.0510 ±0.0010 


16^ 


6.0 


0.82 ± 0.03 


1.0996 ± 0.0006 


16^ 


9.0 


0.91 ±0.04 


1.0460 ± 0.0040 



B. Finite- volume and lattice discretization effects 

Next we analyze discretization and finite-volume ef- 
fects. In this section the discussion will be restricted 
to the propagators for the logarithmic definition. Corre- 
sponding data for ghost-gluon coupling, a^^{q'^), is then 
discussed in the next section. 

When checking lattice discretization artifacts, we fix 
the physical volume such that it approximately equals 
that of all data for different a(/3). In contrast, finite- 
volume effects will be analyzed for a fixed /3 varying the 
lattice size. 

To analyze lattice discretization artifacts for the gluon 
and ghost propagators we compare their renormalized 
dressing functions at /? = 5.8, 6.0 and 6.2. Using the 
respective lattice sizes 16^, 24^ and 32^, the physical vol- 
ume is then roughly V ~ (2.2 fm)*. For the renormaliza- 
tion we chose /i w 3.2 GeV, which we find lies well below 
the momenta where discretization artifacts could affect 
the renormalization. The corresponding data is shown in 
Fig. [TT] suggesting that, with respect to precision of the 
data, lattice discretization artifacts are reasonably small. 

To check finite-volume effects we choose /3 = 6.0 and 
vary the lattice size from 16''', 24^ and 32"^. This has been 
arranged for the data in Fig. [T^] where we compare the 
renormalized gluon and ghost dressing functions. One 



12 




momentum a^q^ 



D" 




12*, p = 9.0 
le'*, p = 9.0 



6 8 10 
momentum a^q^ 



12 14 



16 



FIG. 8: The ratio Cgiuon acc. to (|41|l relating the gluon propagator for the two definitions of the gluon field. Data is for /3 = 6.0 
(left) and /3 = 9.0 (right) on a 12* and 16" lattice. 
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FIG. 9: The bare ghost dressing functions versus the lattice momentum squared. Data is for /3 = 6.0 (left) and /? = 9.0 (right) 
and two lattice sizes. Filled (open) symbols are for the logarithmic (linear) definition. 



clearly sees that finite- volume efi^ects are negligible above 
1 GeV, the momentum where the gluon dressing function 
has its maximum. At 0(1) GeV and below a slight vol- 
ume dependence becomes visible. The overall behavior 
resembles that what has been observed for A*^''"-' in other 
studies. In order to see these effects well below 1 GeV, 
much larger lattices are needed, for example, as those 
studied for the linear definition in [3|. 



C. The running coupling 

The running coupling as(q^) for Yang-Mills theories 
can be defined in various ways. Here we use the coupling 
of ghost-gluon vertex in a particular (minimal) MOM 
scheme (see [Sll, [gOI as well as the more recent papers 
[19I EH). It can be defined in terms of the bare, i.e., 
unrenormalized gluon and ghost dressing functions Zqi 
and Zch as follows: 

ar{e)^^Zo,{a',q')Z'^^{a',q'). (44) 



It is a renormalization-group invariant quantity, i.e., 
shifting the cutoff or transforming the right hand 
side into renormalizcd quantities and changing their sub- 
traction momentum fi within the given MOM scheme 
should not alter a^'^(g^). Therefore, we can compute it 
directly from the bare lattice dressing functions at an ar- 
bitrary large enough cutoff- value a~^(/3), as long as mul- 
tiplicative rcnormalization is ensured and additive lattice 
artifacts are suppressed. In what follows we shall omit 
the superscript MM for simplicity. 

First, we check effects due to the lattice discretiza- 
tion and the finite volume. In order to investigate lattice 
discretization effects we present on the left hand side of 
Fig. [131 the running coupling as(g^) for different lattice 
spacings but fixed physical volume [as above we choose 
again V = (2.2 fm)''']. Apparently, there are some sys- 
tematic lattice discretization effects, suggesting that for 
the given (rather small) (3 additive lattice artifacts are 
small but not negligible. This is in agreement with the 
findings in [l^ . These artifacts should disappear for large 
/3. On the right hand side of Fig. [T31 we show data for 
different physical lattices sizes but fixed lattice spacing 
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FIG. 10: The ratio Cghost acc. to (j42} relating the ghost propagator for the two definitions of the gluon field and the Faddeev- 
Popov operator, respectively, for 12* and 16* lattices generated at /? = 6.0 (left) and /3 — 9.0 (right). 
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FIG. 11: Renormalized gluon (left) and ghost dressing function (right) for the logarithmic definition and various lattice spacings 
a = a(/3). The physical volume is fixed to V = (2.2 fm)*. Data has been renormalized at g = /i ~ 3.2 GcV. 



[again we choose (3 = 6.0]. Based on that figure we have 
to conclude that finite-volume effects seem to be negligi- 
ble for the considered momentum range. 

In Fig. [14] we finally compare the coupling for the log- 
arithmic and the linear definition. We decided to show 
data for various lattice sizes and /? values in a single 
plot to demonstrate that altogether the data for the two 
definitions matches up almost completely. Regrettably, 
there are some small deviations due to the different lat- 
tice spacings and volumes, but this should be expected. 
The almost perfect overlap of the two curves agrees, of 
course, well with what we saw above for the ratios of the 
propagators (see Eq. (|43|)). That is, as{q^) will not dif- 
fer calculated cither for the standard (linear) approach 
(as in [l^ ) or for the logarithmic approach as done here. 
Unfortunately, we cannot say which approach comes with 
the smaller lattice discretization artifacts. This is left for 
a future study. 



VI. COMPARISON WITH NSPT 

We now turn to the NSPT results of [H, [s^, IMI which 
below we will compare to our data. Let us start first with 
some facts on NSPT. 

NSPT is a numerical approach to lattice perturba- 



tion theory that allows to circumvent the difficulties of 
the standard diagrammatic approach and facilitates au- 
tomatized perturbativc calculations. It has its roots 
in stochastic quantization and is based on a modified 
Langevin equation equipped with stochastic gauge fixing 
corresponding to a gauge fixing term (9^A^)^/(2^) at fi- 
nite ^ ^. For our purposes it is actually a hierarchy of 
first-order evolution equations associated with the vari- 
ous parts of the gauge field when expanded in terms of 
the lattice coupling go oc 1/v^: 

U,,,^l + J2r'^'U^%, (45) 
1>1 

From a numerical point of view, these different parts, 
representing first, second, third etc. orders, are separately 
dealt within the code. The maximal addressable order 
of perturbation theory is thus limited by the available 
computing resources (cpu time and memory). 



For a study of an approximate Landau gauge within a Langevin 
approach see Ref. 1631 . 
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FIG. 12: The renormalized gluon (left) and ghost dressing function (right) for the logarithmic definition. The data is for a 
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a — a{j3) but fixed physical volume V = (2.2 fm)''. Right: for different volumes V but fixed lattice spacing a = a(/? = 6.0). 



The Langevin simulation is implemented in an Euler 
scheme with a finite evolution time step. Lattice observ- 
ables, in our case the ghost [1^ and gluon [s^] propaga- 
tors, are evaluated taking the long-time average, order 
by order in a loop expansion in even powers of go. Con- 
tributions from odd powers vanish within the statistical 
errors. As for any Langevin simulation, one then has to 
take the limit to vanishing time step. This is in addition 
to the continuum limit and the limit of infinite volume. 

Regardless of this, NSPT results for a finite lattice vol- 
ume can be confronted directly with standard MC results 
for a given /3, supposed the lattice size and definition of 
the studied observable is the same in both approaches. 

But before such a comparison is possible, the limit 
^ {minimal Landau gauge) has to be taken. This 
is arranged such that a sequence of configurations (sepa- 
rated by 0(50) Langevin time steps) undergoes a Fourier- 
accelerated gauge-fixing procedure, after which the indi- 
vidual gluon fields, A^ii\ each associated with particular 
perturbative order (^q), are transversal within machine 
precision. 

As in standard lattice perturbation theory, an ex- 
panded version of the logarithmic relation between the 



gluon fields and the transporters (compare ([3])) is taken 
into account up to the maximal order of perturbation 
theory addressed in the given case. Correspondingly, the 
gauge functional and the structure of the Faddecv- Popov 
operator are the same as in Eqs. ((37)) and (|38|) . 

The gluon two-point function in n-loop order is then 
defined as a convolution of the bilinears of gluon fields 
(in momentum space) in complementary orders: 

s'''DilHp{k)) - ( E [^;/')(fc)^^<'"+'-'n-fc)] ) ■ 

\ 1 = 1 I u 

(46) 

The Faddeev-Popov operator (explicitly written in 
Eq. (|38p up to fourth order) can be expanded in terms of 
products of various A'^^\ with the term M'"^ collecting 
all terms of order g^. This structure allows to express the 
inverse of the Faddeev-Popov operator also as an expan- 
sion in orders of g^ in a recursive way, without the need 
of explicitly inverting any other than the zeroth order 
term, M(°) = A (the Laplacian). 

A reasonable "convergence" of the NSPT results up 
to few loops (three or four are available now) requires 
a small bare coupling 50 j i-c, a large /3. However, the 
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FIG. 14: Running coupling for various lattice sizes and /?- 
values. Filled symbols are for the logarithmic definition, while 
open symbols are for the linear one. 



bare coupling; is known to be a poor expansion pa- 
rameter [63|. One can speed up convergence by "boost- 
ing", i.e., trading the bare coupling constant by an effec- 
tive "boosted" coupling gl = 5o/-fport(5o) > 5o' where 
Ppert is defined by the average perturbative plaqucttc. 
Its expansion is determined within the Langevin simu- 
lations, along with the propagators. The effect of the 
larger boosted coupling is overcompensated by the rapid 
decay of the expansion coefficients with respect to gl- 

We now compare our data for the logarithmic defini- 
tion to the NSPT results of [H [13, [Ml . 

One should be aware that on small lattices, with 
sizes like 12'* or 16^, and for the larger /3 values, the 
Monte Carlo lattice gauge fields will be in a pseudo- 
deconfinement phase. This can be monitored by the 
"spatially" averaged Polyakov loops, actually in all four 
directions. It is well-known that the agreement of Monte 
Carlo results with standard LPT requires that the Monte 
Carlo simulation is guaranteed to stay in the trivial 
Polyakov sector [3]. For SU{i) this means, that the 
average Polyakov loops in all four directions have to be 
located in the sector of predominantly real values. Let 
us denote averages in this sector by the corresponding 
phases (0,0,0,0), in distinction to results obtained from 
Monte Carlo configurations without taking notice of the 
Polyakov sector. Note that one can easily switch between 
the Polyakov loop sectors by applying global Z{?>) trans- 



formations on all link variables attached to and pointing 
forward, orthogonal to an arbitrary 3D plane. While at 
(3 = 9.0 the MC results for the gluon and ghost dress- 
ing functions clearly depend on the Polyakov loop sector, 
for p ~ 6.0 the Polyakov loop values fluctuate closely 
around the origin of the complex plane, i.e., in the con- 
fined phase. In this case the choice of a sector should 
not influence the behavior of the propagators or dressing 
functions. In what follows, for the detailed comparison all 
Monte Carlo configurations at /3 = 9.0 have been fiipped 
to the real sector (0, 0, 0, 0) if necessary, before the gauge 
fixing has been performed. 

Let us first confront the tree level and the cumulative 
one-loop and two-loop contributions to the gluon dressing 
function with the results from Monte Carlo simulations 
(see Fig. Iisp . The simulations have been performed for 
a 16* lattice at /3 = 6.0 and p = 9.0, the same values as 
for the NSPT results. For the reader's convenience, we 
present the NSPT data the same way the Monte Carlo 
data has been present above. When looking at the data 
in Fig. [15] we see the NSPT results approaches the MC 
data with increasing loop order. As expected though, 
this is less the case for /3 = 6.0, but for /3 = 9.0 the 
NSPT data almost approaches the MC data. 

Similar we see for the bare ghost dressing function in 
Fig. [ini For this the speed of convergence of the NSPT 
results can be assessed from the difference between two- 
loop and three- loop [l!] ■ The three- loop result is already 
very close to the MC result for /3 = 9.0 and at the largest 
available momenta. 

Next we illustrate the effect of "boosting" the pertur- 
bative expansion. For this we use the currently avail- 
able NSPT resuhs (up to four/two loops (L < 12/L = 
16,20,32) for the gluon propagator and up three loops 
[L < 20) for the ghost propagator) and confront them 
with corresponding MC data. This is shown in Fig. [IT] 
where also the bare inverse coupling /3 and its boosted 
value /3boost are given. As expected, boosting moves the 
NSPT data closer to the MC results, but they cannot be 
reached, certainly not at /3 = 6.0. 

Last but not least, the running coupling, as(q^), as cal- 
culated from the NSPT dressing functions, both summed 
up to the available orders, is compared to the Monte 
Carlo results at /? ~ 6.0 and 9.0. The corresponding 
data is shown in Fig. I18[ again for naive and boosted 
perturbation theory. We sec that the running coupling 
from Monte Carlo simulations is approached from below 
up to 7% for P ^ 6.0 and practically approached within 
the present errors for (3 = 9.0. 

Concluding this section, one can say that our MC re- 
sults for the logarithmic approach to Landau gauge prop- 
agators support the validity of the NSPT calculation for 
the gluon as well as for the ghost propagator. In fact, the 
NSPT results do not coincide with the MC data but they 
become closer with increasing order of the perturbative 
expansion. Also, the difference between the NSPT and 
MC data becomes smaller for larger /?, suggesting the 
difference to be related to nonperturbativc effects NSPT 
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FIG. 15: Comparison of MC with NSPT results for the bare gluon dressing function at (cumulative) tree, 1-loop and 2-loop 
level at /3 = 6.0 (left) and ^ = 9.0 (right). The lattice size is 16"*. 



1.45^ 

1.40- 

1.35- 

1.30- 

1.25- 



1-20-^ 



MC-simulation, (0,0,0,0) 
-5- NSPT, 3-loop 

NSPT, 2-loop 
^ NSPT, 1-loop 



1.15-1 



1.10-^ ^^Sb5v„ 
1.05- 



1.00- 



° □ □ 



1.25 



10 12 14 16 



1.20- 



1.15 



o 

N 



1.10 



1.05 



1.00 







MC-simulation, (0,0,0,0) 






NSPT, 3-loop 




-e- 


NSPT, 2-loop 


w 




NSPT, 1-loop 


w *■ 




























A 




" 


^ W * ** i 

X V ^ ^ 



momentum a^q^ 



6 8 10 
momentum a^q^ 



12 14 16 



FIG. 16: Comparison of MC with NSPT results for the bare ghost dressing function at (cumulative) 1-loop, 2-loop and 3-loop 
level at /3 = 6.0 (left) and /? = 9.0 (right). The lattice size is 16*. 



results cannot provide. 

VII. CONCLUSIONS 

In this paper we have studied an alternative approach 
to compute the SU{3) Landau gauge gluon and ghost 
propagators on the lattice. This approach uses a logarith- 
mic ansatz [Eq. (|3])] for the definition of the lattice gluon 
fields from a given gauge field configuration. It is thus 



best suited to compare lattice MC data for these propaga- 
tors with results from NSPT. We have started the task by 
first exploring some options for an efficient algorithm that 
fixes gauge field configurations to Landau gauge for the 
logarithmic case. We find a multigrid-accelerated gra- 
dient method with a preconditioning step of simulated 
annealing and subsequent overrelaxation, applied to the 
linearly defined gluon field, to be a good choice. The 
method is also easy parallelizable. 

With this algorithm at our disposal we have com- 
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pared the bare lattice propagators for the logarithmic 
and the linear definition of gluon fields. As expected, we 
see them to differ by multiplicative factors that depend 
on /3. Those factors are such that they perfectly can- 
cel when one considers Q;s(g^). Apart from some small 
lattice discretization artifacts, data for the running cou- 
pling matches up almost completely, as it should for a 
rcnormalization-group invariant object. It is thus also 



an ideal quantity to assess discretization artifacts. 

For the logarithmic definition we have checked Gribov 
copy, finite-size and lattice discretization effects, and find 
them to be small for momenta q> \ GcV. 

Finally we have compared our MC data for the log- 
arithmic definition with results from NSPT. These are 
available up to four loops for the gluon propagator, and 
up to three loops for the ghost propagator. Wc find a rea- 
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sonable convergence at large momentum. Note that for 
this it is important that during the MC process the gauge 
field configurations are kept in the correct (real-valued) 
Polyakov loop sector. For large f3, these may easily pass 
a pseudo-deconfinement phase transition. 

Our results altogether support universality with re- 
spect to the two lattice realizations of SU{3) Landau 
gauge theory studied herein. In as far this universality 
persist in the low-momentum region remains to be seen. 

We emphasize that the universality of different lat- 
tice definitions we have observed in this paper assumes a 
unique Landau gauge fixing based on the (global) min- 
imization of a corresponding gauge functional. For an 
alternative view of dealing with the Gribov copy prob- 
lem see Refs. HEil. 



Multigrid of lattice size h' 



calculate Eq. (A3) with lattice spacing h' 
if lattice h' is coarsest then 



calculate projection: 



h' 



Eq. (A5) 



> h- 



iteratively start p times 



Multigrid of lattice size h 



1 1 • 1 • , Eq. (A4) , , 

calculate mterpolation: n > h ~ 

calculate Eq. (A3) with lattice spacing /i'- 
output: w'' 



Acknowledgments 



FIG. 19: Flow chart illustrating our multigrid algorithm. 



The authors are grateful for the generous support by 
the HLRN Berlin/Hannover with CPU time on their 
massively parallel supercomputing system. E.-M.L and 
M.M.-P. acknowledge financial support by the DFG via 
Mu932/6-l. A. St. acknowledges support from the Eu- 
ropean Union through the FP7-PEOPLE-2009-RG pro- 
gram and by the SFB/Tr 55. 



Appendix A: Multigrid Fourier-accelerated 
gauge-fixing 

For our implementation of a parallel version of the 
Fourier-accelerated gauge fixing we follow Goodman and 
Sokal [Sll, [13] using the representation of the Fourier 
transformation of l/<z^ in position space: 



F' 



ymax 



F 



(^9/4'°s))(, 



(Al) 



This leads us to a simple inversion of the Laplacian, 
i.e., to solving the 4D Poisson equation 



(a«)(.t) = (x) 



(A2) 



using for example the local Jacobi method. To avoid 
critical slowing down inherent to this method we use the 
multigrid algorithm by solving ()A2|) successively on the 
original (fine) lattice and on several coarser lattices. In 
order to use our multigrid algorithm for various lattice 
sizes we implement the multigrid with a symmetric lattice 
decomposition. 

On each sublattice one has to solve a system of linear 
equations 



(A3) 



where the superscript h describes the respective lattice 
spacing h = a, 2a, 4a, 8a, ... Comparing with \K2\ we can 
see that = A. For switching between the lattices one 
defines interpolation matrices / as 



j^h' _ ^ Ah I T iT 



h'.hA'^ [Ih',h] 



= h'Mv" , = Ih',hf 



(A4) 



and projection matrices P 



Ph.h'V 



h' 



f' = Ph,h'.f 



h' 



(A5) 



h' always denotes the lattice with the finer spacing while 
h denotes the coarser lattice. The matrix structure (in 
terms of the lattice sites) of these equations is left im- 
plicit. The projection matrices are the transposed inter- 
polation matrices (with respect to the indices pointing to 
lattice sites) 



Phji' — [Ih'.h]^ 



(A6) 



The matrices were chosen in the way that the operator 
got the same structure on all sublattices (i.e. A^ = 
A'' ). In practice, the multigrid algorithm is realized by 
jumping between the finest and various coarser lattices. 
This is summarized in the flow chart in Fig. 1191 

Solving Eq. (|A3p before the projection is called pre- 
smoothing, post-smoothing after the interpolation, re- 
spectively [65|- To solve Eq. (jA3p we use the Jacobi 
method with a fixed number of 20 iterations. Hence we 
did not solve equation Eq. (jA3[) at high numerical ac- 
curacy on the sublattices. Nevertheless, the accuracy of 
this calculation seems to have only a minor influence on 
the total number of iterations needed to fix the Landau 
gauge. For the parameter p we chose the value p = 2, 
i.e., the so called VF-cycle. 
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